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I.  INTRODUCTION 


The  theory  and  technique  for  determining  the  temperatures  of  gases  from 
their  Fourier  transform  infrared  spectra,  in  emission  and  absorption,  has  been 
well  documented  in  the  literature.  Frequently,  temperature  determination 
involves  fitting  peak  absorbances  of  rotational  fine  structure  within  a 
vibrational  transition  to  a  Boltzmann-type  population  expression,-  since  the 
relative  absorbances  are  strongly  temperature  dependent.  One  of  the 
difficulties  associated  with  calculating  the  temperature  of  a  gas  from  the 
peak  absorbances  of  ro-vibrational  transitions  is  instrumental  distortion  of 
peak  absorbance  values.  This  distortion  results  from  the  convolution  of  the 
source  intensity  (modulated  by  the  sample)  and  the  instrument  lineshape 
function  (ILS),  and  is  most  important  in  spectra  obtained  at  moderate  to  low 
resolution  (here  taken  to  correspond  to  optical  retardations  of  less  than 
~1.0  cm).  The  method  for  correcting  the  error  in  observed  peak  absorbance  due 
to  this  convolution  as  applied  to  FTIR  spectroscopy  was  first  described  by 
Griffiths. 


The  ana] ytical  treatment  of  most  physically  encountered  convolutions  is 
difficult.  To  our  knowledge,  no  analytical  solution  to  the  resulting 
"convolution  integral"  for  transmission  spectra  obtained  on  a  FTIR 
spectrometer  has  been  reported.  Ue  have  recently  investigated  and  quantified 
the  relationship  between  observed  and  undistorted  peak  absorbances  measured  by 
FTIR  spectroscopy.  This  relationship  was  first  quantified  by  Griffiths  and 
Anderson.  Our  goal  is  to  develop  a  reasonably  accurate  (within  ~5%) 
graphical  method  to  obtain  the  temperature  of  a  gas  sample  at  elevated 
temperatures  (>500  K)  from  absorption  spectra  in  the  infrared  using  resources 
typically  available  in  an  analytical  laboratory.  We  assume  this  to  be  a 
Fourier  transform  IR  spectrometer  capable  of  1  cm-*  resolution  and  a 
microcomputer  of  virtually  any  type.  The  graphical  method  is  not  intended  to 
supplant  the  more  accurate  methods  already  in  the  literature , ^ but  rather  to 
provide  a  fast  ano  convenient  complementary  view  of  the  data.  To  date,  we  are 
unaware  of  any  successful  temperature  estimations  from  peak  absorbances  of  ro- 
vibrational  transitions  which  have  employed  graphical  methods  over  a 
temperature  range  of  several  hundred  degrees. 


II.  BACKGROUND 


For  a  Fourier  transform  spectrometer  employing  triangular  apodization  of 
the  interferogram,  the  observed  minimum  transmittance,  T0,  of  an  absorption 
feature  with  a  Lorentzian  profile  centered  at  frequency  vu,  may  be  given  by:"1 

Tn  =  _l  exp{-2.303Apma^/((v-vn)2+a^'}Dsin2[piD(v-vI3)]/ 

[piD(v-v]3)  ] 2  dv 


where  A„„  is  the  true  weak  absorbance  of  line  d  at  freauencv  v.,.  a  is  th» 
linewidth,  D  is  the  maximum  optical  retardation  of  the  spec' rometer,  and  pi  is 
3.1416.  Here,  ra  is  the  line  index  number,  where, m=-J"  for  P  branch 
transitions  and  n=J"+l  for  R  branch  transitions.0  This  equation  describes  the 
minimum  transmittance  of  an  isolated  pressure  broadened  line  where  the 
absorption  feature  is  unobscured  by  an  neighboring  lines,  equation  (1)  is  a 
function  of  temperature  through  the  dependence  of  t/ie  linewidth,  an,  on 


1 


temperature.  The  variation  of  a^  with  temperature  and  pressure  may  be 
calculated  by  the  approximation: 

am(T,P)  =  an0(P/P0)  (T0/T)n  (2) 

where  Pq  and  Tg  are  a  pressure  and  temperature  for  which  the  linewidth  a,3g  is 
known,  P  and  T  are  the  experimental  pressure  and  temperature,  and  n  is  the 
empirically  determined  temperature  broadening  exponent.  The  coefficient  n  has 
been  shown8  to  be  temperature  and  line  dependent.  For  the  sake  of 
simplification,  n  is  usually  set  to  0.75,  although  there  is  evidence8  that  a 
value  of  0.66  may  be  more  accurate  at  room  temperature. 

9 

True  peak  absorbance  may  be  fit  to  temperature,  T,  using: 

Apm  a  |ra|F(m)vrlexp{-E(m)/kT} /am  (3) 

where  au  varies  as  in  Eq.  (2),  E(m)  is  the  rotational  energy  of  the  lower 
state  in  the  transition.  F(m)  is  the  Hermann-Wallis  factor  for  interaction  of 
rotation  and  vibration,  and  |n|  is  the  absolute  value  of  the  line  index 
number.  The  linewidth  appears  in  the  denominator  of  Eq.  (3)  to  take  into 
account  the  dependence  of  the  peak  height  on  ygak  width  for  a  spectral 
absorption  feature  with  a  Lorentzian  profile. 

By  substituting  Eq.  (2)  into  Eq.  (1)  and  obtaining  an  analytical  solution 
to  the  integral,  an  expression  relating  observed  minimum  transmittance  (and 
hence  observed  peak  absorbance)  to  true  peak  absorbance  may  be  obtained.  If 
we  denote  this  expression  as: 

V  =  £(Van.D>  <4) 

and  combine  this  with  Eq.  (3)  we  obtain: 

f(Tra,au,D)  a  | m| F(m)vnexp{-E(n)/kT} I aQ  (5) 

The  first  step  in  linearizing  Eq.  (5)  involves  finding  an  approximation 
to  Eq.  (1).  Griffiths  gives  the  expression:  ’ 


A  a  A„K c  W  (6) 

pm  obs 

where  R,  the  resolution  parameter,  is  equal  to  l/(2Dan)  and  x  and  y  are  set 
equal  to  2.  If  this  expression  is  substituted  into  Eq .  (5),  and  the 
temperature  broadening  exponent  n  is  assumed  to  be  invariant  with  m,  one 
obtains: 

AobsX/ f  |n|F(n)vm)}  a  exp{-E(n)/kT}  (7) 

When  ambient  temperature  spectra  are  collected  (for  which  the  a,3  for  many 
diatomics  are  accurately  known)1  and  the  resolution  parameter  for  the  ro- 
vibrational  lines  being  used  in  the  calculation  is  such  that  there  is  a  near 
linear  relationship8  between  log(A  )  and  log(A0^s),  a  graph  of  E(m)/k  vs. 
ln(A  bsX/{  (2Darlo)y_1(  |tn|F(m)vi;i)}  ]  produces  a  straight  line  with  a  slope  (1/T) 
the  inverse  of  which  usually  falls  within  a  percent  of  the  true  gas 
temperature As  the  temperature  of  the  gas  being  investigated  is  increased, 
the  error  between  the  calculated  and  actual  temperature  increases.  This 
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increase  in  error  is  caused  by  the  change  in  an  with  T,  which  affects  the 
resolution  parameter  R,  and  nay  shift  the  relationship  between  log(Apn)  and 
log(AQks)  into  a  nonlinear  region.  Also,  as  stated  earlier,  the  temperature 
broadening  exponent,  n,  from  Eq.  (2)  changes  with  temperature  and,  for  a  given 
temperature,  nay  vary  with  line  index  number. 

Til.  CALCULATIONS 


Ue  have  fir  input  parameters  and  the  result  of  the  numerical  integration 
of  Eq.  (1)  to  an  equation  of  the  form: 


^■pm  “  k^obs  ^ 


(8) 


where  k,  x,  and  y  are  constants  and  R  is  the  spectral  resolution  parameter. 

The  integral  in  Eq.  (1)  was  evaluated  numerically  for  all  combinations  of 
optical  retardzations  of  1,  2,  4,  8,  and  10  cm,  true  peak  absorbances  of  0.1 
to  10  (in  intervals  of  0.1)  with  absorbing  lines  for  which  an  ranged  from  0.01 
to  0.1  cm-i  (in  intervals  of  0.01  cm-*).  A  weighted  line  -  ieast  squares 
routine  was  used  to  fit  the  data  from  portions  of  the  numerical  integration  of 
Eq.  (1)  to  a  linearized  form  of  Eq.  (8). 

The  integration  routine  assumed  that  vn  was  located  at  2000  cm  and 
elevated  Eq.  (1)  using  Simpson' s1^  rule  for  50  cm-1  on  either  side  of  vn*  The 
interval  to  be  numerically  integrated  was  divided  into  10,000  segments.  The 
choice  of  value  for  the  center  of  absorption,  vm,  is  arbitrary,  as  may  be 
verified  by  substituting  v-vQ  in  Eq.  (1).  To  determine  if  numerical 
integration  gave  values  dependent  on  the  value  of  vn  ranging  from  500  to 
3500  cm  1  showed  the  value  of  the  integral  to  be  independent  of  the  value  for 
this  parameter. 


The  application  of  our  numerical  integration  routine  to  the  function 
sin  .x/x^  (which  is  similar  in  behavior  to  Eq.  (1)  yet  has  an  analytical 
solution)  using  parameters  similar  to  those  in  the  numerical  integration  of 
(1)  gave  results  within  0.001%  of  the  analytical  solution.  Ue  believe  our 
numerical  integration  of  Eq.  (1)  gives  similar  accuracy. 

For  all  ranges  of  Apn,  values  used  in  the  integration  of  Eq.  (1),  the 
average  deviation  between  the  numerically  integrated  value  -log(Tn)  and  the 
fitted  approximation  presented  here  was  less  than  that  of  Eq.  (6)  with  x  and  y 
equal  to  2.  In  fairness,  the  authors  of  Reference  5  presented  their 
approximation  without  derivation  and  in  no  way  implied  it  to  be  valid  over  any 
range  other  than  that  specified  in  the  paper.  In  addition,  for  the  range  over 
which  the  approximation  of  Reference  5  was  used,  that  approximation  and  the 
fitted  approximation  developed  here  are  in  good  agreement. 

For  spectra  of  gases  below  800  K,  the  integral  in  Eq.  (1)  was  evaluated 
for  true  peak  absorbances,  Apn,  ranging  from  1  to  5  absorbance  units,  an 
optical  retardation  of  1  cm,  and  an  values  ranging  from  0.01  to  0.1  cm”  . 

These  pa i duelers  are  typical  for  CO  at  pressures  of  several  hundred  torr. 

The  best  fit  values  obtained  for  In(k),  x,  and  y  in  Eq.  (8)  were  0.122  (.014), 
1.818  (.014),  and  1.748  (.015),  respectively,  where  the  numbers  in  parentheses 
are  equal  to  one  standard  deviation.  For  the  same  parameters  mentioned  above, 
but  for  peak  absorbances  ranging  from  0.1  to  1.0  absorbance  units,  the  best 
fit  values  of  ln(k),  x,  and  y  were  respectively,  0.211  (.00004),  1.287 
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(.00002),  and  1.208  (.00002).  The  spectroscopic  constants  used  to  calculate 
E(n)  and  F(n)  nay  be  found  in  Reference  13.  Values  of  vn  used  in  the 
calculation  are  those  reported  by  the  spectrometer .  Values  of  an  were 
obtained  from  Reference  5.  For  this  work,  it  was  assumed  that  ail  gases  had 
the  sane  broadening  effect  on  the  am  for  CO  as  did  CO.  The  additional  error11 
probably  is  insignificant  given  the  assumption  that  n  is  constant. 


At  temperatures  below  600  K,  20  lines  from  the  R  branch  (because  it  is 
more  intense  than  the  P  branch  and  for  reasons  given  below)  in  the  fundamental 
absorption  band  of  CO  were  used  in  the  calculation.  To  calculate  the 
temperature,  Eq.  (7)  was  first  linearized,  and  a  v.  ghted  linear  least 
squares1^  routine  used  to  determine  1/T. 


To  decide  which  lines  will  be  used  to  calculate  temperatures  from  FTIR 
spectra,  it  is  necessary  to  establish  a  criteria  for  line  acceptability. 

Lines  which  are  significantly  overlapped  by  absorbances  due  to  hot  band 
transitions  or  to  absorbances  by  isotopically  substituted  honologues  should  be 
excluded  from  the  computation.  Figure  1  shows  a  portion  of  the  CO  spectrum  in 
a  cell  held  at  866  K.  Visible  in  this  figure  are  absorbances  due  to  1-1C0  and 
those  due  to  the  1-2  transition  in  1  CO,  as  well  as  the  much  more  intense 
fundamental  transition.  The  error  caused  by  hot  band  overlap  is  very 
temperature  sensitive  and  becomes  most  important  at  high  temperature.  Figures 
2  and  3  show  the  line  center  separation  in  wavenumbers  between  fundamental  and 
hot  (v=l-2)  band  transitions  in  ordinary  CO  between  fundamental  transitions  in 
normal  and  3C  substituted  CO,  respectively.  These  figures  were  generated  by 
calculating  rotational  energy  levels  from  spectroscopic  constants1-1  for  the 
appropriate  species/vibrational  level,  calculating  transition  frequencies 
using  conventional  selection  rules,1-3  ordering  these  transitions  with  respect 
to  the  energy  of  transition,  and  then  taking  the  difference  in  energy  between 
each  successive  level.  This  energy  deference  is  plotted  versus  the  energy  of 
the  transition  with  the  lower  energy.  For  each  species  considered,  J"  values 
from  0  or  1  to  30  were  used.  At  a  glance,  each  of  these  figures  provides  a 
good  idea  of  where  regions  of  minimum  line  overlap  exist  for  the  fundamental 
absorption  of  CO.  Because  these  figures  doe  not  include  intensities  of 
individual  transitions,  they  cannot  indicate  where  a  weak  fundamental 
absorption  line  used  in  the  calculation  is  overlapped  by  a  signf iciantly 
strong  ^CO  line,  for  example,  or  vice  versa.  For  CO,  typical  ranges  of  lines 
to  use  in  the  calculation  are  the  20  lines  in  the  regions  2059-2139  cm-1  (P- 
branch  J"=l-20)  and/or  2147-2212  cm  1  (P.-branch  J"0-19).  An  examination  of 
Figure  2  shows  that  these  two  regions  of  the  fundamental  absorption  in  CO 
significantly  overlap  with  the  v=l-2  transitions  in  CO,  with  the  P  branch 
being  slightly  less  obscured.  However,  since  the  1-2  envelope  is  shifted  by 
about  25  cm  to  lower  energy  than  the  0-1  envelope  (not  apparent  from  the 
figure),  tne  P.  branch  of  the  0-1  transition  is  relatively  unobscured  by  those 
1-2  transitions  which  have  appreciable  intensity  at  normal  temperatures. 


Figure  3  shows  that  both  branches  in  the  fundamental  absorption  of  CO 
have  several  lines  with  separation  between  line  centers  for  overlap  with  1-1C0 
of  less  than  0.5  cm  .  However,  because  the  band  centers  differ  between  CO 
and  CO  by  about  45  cm  ,  the  R  branch  of  CO  is  relatively  unobscured  by 
overlap  from  the  R  branch  of  ^CO  (since  the  R  branch  of  ^C0  overlaps  the  P 
branch  of  CO).  Therefore,  Figures  2  and  3  suggest  using  the  R-branch  of  CO 
foe  temperature  determinations  at  temperatures  where  hot  band  transitions  are 
not  important.  At  higher  temperatures,  Figure  3  suggests  a  switch  to  the 
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Figure  3.  Line  Center  Separation  Between  the  v=0-l  Transition  in  CO  and  the 

v=0-l  Transition  in  L'}C0  (see  text) 


P-branch  of  CO  for  temperature  determinations  when  1-2  transitions  become  more 
troublesome  than  JC0  transitions. 

IV.  EXPERIMENTAL 

All  temperatures  reported  here  were  obtained  from  absorption  spectra. 
Since  the  heated  cell/burner  was  located  in  the  sample  compartment  of  an 
otherwise  unmodified  FTIR,  infrared  emission  was  unmodulated,  and  no  attempt 
was  made  to  collect  emission  spectra.  No  correction  was  made  to  the  spectra 
to  account  for  emission.  Experiments  were  performed  on  a  Mattson  Sirius  100 
FTIR  spectrometer  using  a  HgCdTe  liquid  nitrogen  cooled  detector.  All  spectra 
were  obtained  with  100  co-added  scans  obtained  at  1  cm  resolution. 

Interf erograms  were  triangularly  apodized  and  zero-filled  to  the  limit  of  our 
computer  memory  (131072  points)  prior  to  transformation. 

Temperature  controlled  stainless  steel  cells  were  ec  lipped  with  NaCl 
windows.  The  pathlength  for  absorption  in  each  cell  was  pproximately 
12  cm.  For  the  temperature  range  300  K  to  A 7 3  K,  the  NaCl  windows  were 
attached  to  the  cell  body  using  metal  flanges  and  fiber  gaskets.  The  leak 
rate  of  this  cell  at  all  temperatures  used  was  always  less  than  500  mtorr/ 
hr.  Pressures  used  in  these  low  temperature  experiments  (<473  K)  ranged  from 
115  to  197  torr.  For  the  temperature  range  from  473  K  to  773  K.  the  NaCl 
windows  were  attached  directly  to  the  cell  body  by  pressure  supplied  by 
flanges.  To  prevent  heat  loss,  the  high  temperature  cell  was  encased  in 
asbestos  which  extended  approximately  1  inch  past  the  windows  on  either  side 
of  the  cell.  No  other  precautions  were  taken  to  avoid  temperature  gradients 
along  the  optical  pathlength.  The  leak  rate  of  this  cell  wa  such  that  all 
high  temperature  experiments  were  performed  at  760  torr  total  pressure  using  a 
4:1  mixture  of  air:C0.  A  check  of  the  instrument  after  filling  this  cell, 
allowing  it  to  stand  for  a  typical  time  of  one  experiment  (~10  minutes)  and 
then  evacuating  the  cell  gave  not  CO  absorption  due  to  CO  in  the  optical 
compartment  of  the  spectrometer. 


Temperature  was  controlled  by  an  Omega  Model  2010  temperature  controller 
using  heating  cartridges  placed  in  cavities  in  the  cell  body  and  monitored  by 
copper-constantin  thermocouples  placed  near  the  inner  wall  of  the  cell  and 
placed  within  the  gas  volume  but  out  of  the  optical  path.  A  check  of  the  cell 
temperature  by  an  Omega  Model  2168A  digital  thermometer  inserted  into  the  body 
of  the  cell  gave  temperature  readings  within  5  K  of  those  reported  by  the 
controller  thermocouple.  All  temperatures  reported  as  actual  are  the 
temperatures  of  the  cell  body.  After  filling  the  cell,  the  gas  was  allowed  to 
stand  for  three  to  five  minutes  to  reach  the  temperature  of  the  cell.  Waiting 
for  longer  periods  of  tine  gave  no  significant  change  in  calculated 
temperature.  No  correction  was  made  for  temperature  variations  in  the 
vicinity  of  the  cell  windows.  At  temperatures  greater  than  550  K,  some 
conversion  of  CO  to  CO2  was  observed  for  the  cell  containing  CO  plus  air. 


Tlie  experimental  setup  for  the  low  pressure  flame  experiments  is  shown  ir 
Figure  4.  The  optical  pathlength  within  the  burner  chamber  was  approximately 
6  cm.  The  burner  head  was  constructed  entirely  of  Pyrex,  with  the  burner 
housing  constructed  of  brass  and  aluminum.  Gas  flow  to  the  mixing  section 
below  the  burner  head  was  regulated  by  calibrated  floating-ball  type 
flowmeters.  Pressure  was  measured  using  an  MK5  Instruments  Baratron 
capacitance  manometer  and  regulated  by  adjusting  a  valve  downstream  from  the 
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Capacitance  Manometer 


Figure  4.  Experimental  Setup  for  Burner  Flame  Experiments 
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flowmeters.  Windows  were  of  LiF .  A  low  pressure  (<100  torr  operation),  the 
flane  is  lifted  off  the  burner  surface  by  2-4  nm,  preventing  damage  to  the 
Pyrex  burner  frit.  Data  was  collected  at  the  maximum  iris  setting  of  the 
spectrometer.  Although  maximizing  sample  compartment  aperture  increases  the 
likelihood  of  resolution  errors  due  to  beam  divergence,  the  bean  waist  at  all 
times  was  kept  within  the  limits  built  into  the  spectrometer.  No  correction 
v/as  applied  to  reflect  any  non-collimated  beam  related  errors. 

V.  RESULTS  AND  DISCUSSION 


The  results  of  temperature  calculations  using  the  temperature  controlled 
gas  cells  for  the  range  300  K  to  773  K  are  shown  in  Table  1.  Values  of  x  and 
y  used  in  calculations  over  this  temperature  range  were  1.818  and  1.748, 
respectively.  At  temperatures  above  600  K,  the  P  branch  of  the  ro-vibrational 
envelope  of  CO  was  used  in  the  calculations  (see  discussion  above).  Figure  5 
shows  a  plot  of  E(n)/k  vs.  ln(A0^sx/ { (2DanjQ)^  Ju|f(u)vu})  for  samples  of  CO 
held  in  the  temperature  controlled  cells  over  the  temperature  range  302  K.  to 
773  K.  The  slope  of  the  best  straight  line  through  each  set  of  data  obtained 
by  the  weighted  linear  least  squares  routine1  is  equal  to  1/T.  In  all  cases, 
calculations  yield  gas  temperatures  within  10%  of  the  cell  temperature. 
Neglecting  the  one  calculated  point  at  620  K,  calculated  temperatures  are 
within  5%  below  700  K.  The  best  results  were  obtained  when  the  temperature  of 
the  gas  was  near  400  K.  This  is  in  contrast  to  earlier  work*1  in  which  the 
best  results  were  found  at  room  temperature.  The  reason  for  this  difference 
is  that  our  fit  to  Eq.  (1)  is  more  accurate  for  small  values  of  Apn  (<3)  than 
for  cases  where  x  and  y  in  Eq.  (6)  are  equal  to  2.  As  average  peak 
absorbances  of  individual  lines  decrease  with  increasing  temperature  (total 
number  of  lines  increases),  an  approximation  to  Eq.  (1)  which  emphasizes  small 
values  of  Apn  yields  a  temperature  that  is  closer  to  actual  than  when  x  and  y 
in  Eq.  (6)  are  equal  to  2. 


The  curvature  in  the  far  right-hand  side  of  the  plots  in  Figure  5  in  the 
region  corresponding  to  low  values  of  J"  is  analogous  to  that  seen  in 
Reference  5.  This  arises  from  the  fact  that  CO  fundamental  absorbances  with 
low  values  of  J"  have  the  larges  afl  and  smallest  value  of  the  resolution 
parameter,  R.  Spectra  of  lines  with  small  resolution  parameters  and  small 
absorbances  have  been  shown^’^  not  to  adhere  to  a  linear  relationship  between 
log(A  )  and  log(AQj3S)  when  obtained  using  triangular  or  sine  instrument  line 
shapeK(IL5)  functions. 


Figure  6  shows  a  portion  of  the  absorbance  spectrum  of  a  rich  low 
pressure  (64  torr)  preraixed  methane/nitrous  oxide  burner  flame  used  in  the 
calculation  of  rotational  temperatures  of  CO.  The  CO  hot  band  (v=l-2)  P 
branch  transitions  were  chosen  for  the  analysis  because  the  hot  band 
transitions  are  shifted  to  lower  energy  than  those  transitionsof  the 
fundamental.  Also,  interference  from  CO2  bands  around  2349  cm  and  from  1^0 
bands  around  2223  cm-1  makes  finding  suitable  lines  in  the  fundamental  of  CO 
difficult.  Hot  band  transitions  were  chosen  to  discriminate  against  any  cold 
CO  near  the  inner  surfaces  of  the  burner  chamber  windows.  Peaks  were 
identified  by  calculating  the  peak  locations  for  the  v=l-2  ro-vibrational 
envelope  of  CO  from  spectroscopic  constants.  i  All  calculations  agreed  within 
a  few  hundredths  of  a  wavenumber  with  the  peak  locations  reported  by  the 
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Hattson  peak  picking  software.  Table  n  shows  calculated  CO  rotational 
temperatures  as  a  function  of  trie  distance  of  the  center  of  the  IR  bean  above 
the  burner  frit  surface. 


Table  1.  Calculated  and  Actual  CO  Temperatures  for  the  Temperature 
Controlled  Gas  Cells  Used  in  These  Experiments  Over  the  Range 

302  to  773  K 


Actual 

Calculated 

Standard 

Total 

Actual-Observed 

mperature 

Temperature 

Deviation 

Pressure 

Temperature 

Kelvin) 

(Kelvin) 

(Kelvin) 

(torr) 

(Kelvin) 

302 

299 

4 

166.8 

3 

302 

297 

4 

164.7 

5 

302 

296 

4 

166.8 

6 

313 

306 

5 

760 

7 

322 

313 

3 

196.7 

9 

323 

319 

5 

163 

4 

325 

323 

6 

141.8 

2 

332 

324 

4 

184.7 

8 

343 

335 

4 

195.1 

8 

343 

337 

4 

189.6 

6 

344 

340 

6 

152.1 

4 

363 

364 

7 

127.3 

-1 

364 

360 

5 

190.3 

4 

364 

360 

4 

193.5 

4 

373 

369 

5 

194 

4 

373 

369 

5 

196.2 

4 

374 

373 

6 

164 

1 

413 

415 

9 

140.9 

_2 

413 

416 

9 

141.4 

-3 

432 

432 

8 

164.8 

0 

433 

436 

11 

115.1 

-3 

438 

437 

8 

184.6 

1 

473 

477 

3 

760 

-4 

473 

478 

5 

760 

-5 

576 

602 

9 

760 

-26 

592 

594 

16 

760 

-2 

674 

620 

61 

760 

54 

674 

672 

13 

760 

0 

t- 

674 

673 

40 

760 

1 

773 

75S 

23 

760 

15 

773 

813 

48 
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-40 
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817 

35 
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-44 
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Figure  6.  Spectrum  of  a  Rich  CH^/f^O  Flame  at  64  Torr.  Peaks  marked 
with  *  are  due  to  the  P  branch  of  the  v=l-2  transition  in  CO. 


Table  2.  Calculated  CO  Rotational  Temperatures  as  a  Function  of  the 
Distance  Between  the  Center  of  the  IR  Beam  and  the  Burner  Frit  Surface 


Burner  Surface 
to  Bean  Center 
(cm) 


Observed 

Temperature 

(Kelvin) 


Standard 

Deviation 

(Kelvin) 
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Figure  7  shows  a  plot  of  E(n)/k  vs.  ln(A0^sx/{ (2Dan  q)^-* |n|F(n)vm) )  for 
the  burner  flame  experiments.  Since  the  observed  absorbances  are  very  weak 
(maximum  less  than  0.12),  wa  used  the  approximation  to  true  peak  absorbance 
from  observed  peak  absorbance  fnv  the  range  of  true  peak  absorbance  from  0.1 
to  1.0  (x=1.287,  y=1.208).  The  flame  temperatures  reported  here  are  lower 
than  adiabatic  fcr  the  nethane/nitrous  oxide  flame  system.  Ue  believe  this 
is  because  the  small  burner  chamber  acted  as  a  heat  sink  for  the  flame.  Mo 
thermocouple  measurements  were  made  of  the  flame  temperature  within  the  burner 
chamber.  Also,  because  we  were  unable  to  find  published  data  for  the  an  of 
the  hot  band  spectral,  we  used  the  values  for  the  fundamental.  Finally^  we 
made  no  provision  for  broadening  of  the  CO  lines  by  gases  other  than  CO,  so  in 
effect,  we  were  treating  our  data  as  if  it  were  obtained  in  a  chamber  of  pure 
CO.  Our  reasons  for  this  stemmed  from  the  many  different  species  present  in 
the  flame  and  our  inability  at  this  tine  to  quantify  many  of  these  species, 
particularly  nitrogen. 


VI.  CONCLUSION 

Calculations  of  rotational  temperatures  of  CO  at  high  temperatures  and  of 
CO  participating  in  combustion  reactions  by  use  of  simple  graphical  methods 
applied  to  FTIR  spectra  have  been  shown  to  be  accurate  within  10%  (and  usually 
much  less)  of  actual  temperature  up  to  800  K.  Graphical  analysis  of  the  ro- 
vibrational  structure  of  the  fundamental  absorption  of  CO  gives  results  which 
indicate  :hat  the  simple  method  employed  here  has  practical  applications, 
especially  when  an  accurate  approximation  to  true  peak  absorbances  is  used. 
Results  from  spectra  of  flames  are  encouraging  but  more  work  is  needed.  Ue 
are  presently  engaged  in  pursuing  this  area  of  research. 
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ATTN:  M.C.  Lin 

J.  McDonald 
E.  Oran 
J.  Shnur 

RJ.  Doyle,  Code  6110 
Washington,  DC  20375 

1  Commanding  Officer 
Naval  Underwater  Systems 

Center  Weapons  Dept. 

ATTN:  R.S.  Lazar/Code  36301 

Newport,  RI  02840 

2  Commander 

Naval  Weapons  Center 
ATTN:  T.  Boggs,  Code  388 

T.  Parr,  Code  3895 
China  Lake,  CA  93555-6001 

1  Superintendent 

Naval  Postgraduate  School 
Dept,  of  Aeronautics 
ATTN:  D.W.  Netzer 

Monterey,  CA  93940 

3  AL/LSCF 

ATTN:  R.  Corley 

R.  Geisler 
J.  Levine 

Edwards  AFB,  CA  93523-5000 

1  AL/MKPB 

ATTN:  B.  Goshgarian 

Edwards  AFB,  CA  93523-5000 

1  AFOSR 

ATTN:  J.M.  Tishkoff 

Bolling  Air  Force  Base 
Washington,  DC  20332 

1  OSD/SDIO/UST 
ATTN:  L.  Cavcny 

Pentagon 

Washington,  DC  20301-7100 

1  Commandant 
USAFAS 

ATTN:  ATSF-TSM-CN 

Fort  Sill,  OK  73503-5600 

1  FJ.  Seiler 

ATTN:  S.A.  Shackleford 

USAF  Academy,  CO  80840-6528 
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1  University  of  Dayton  Research  Institute 
ATTN:  D.  Campbell 

AL/PAP 

Edwards  AFB,  CA  93523 
1  NASA 

Langley  Research  Center 
Langley  Station 

ATTN:  G.B.  Northam/MS  168 

Hampton,  VA  23365 

4  National  Bureau  of  Standards 
ATTN:  J.  Hastie 

M.  Jacox 
T.  Kashiwagi 
H.  Scmcrjian 

US  Department  of  Commerce 
Washington,  DC  20234 

1  Aerojet  Solid  Propulsion  Co. 

ATTN:  P.  Micheli 

Sacramento,  GA  95813 

1  Applied  Combustion  Technology,  Inc. 

ATTN:  A.M.  Varney 

P.O.  Box  17885 
Orlando,  FL  32860 

2  Applied  Mechanics  Reviews 
The  American  Society  of 

Mechanical  Engineers 
ATTN:  R.E.  White 

A.B.  Wenzel 
345  E.  47th  Street 
New  York,  NY  10017 

1  Atlantic  Research  Corp. 

ATTN:  M.K.  King 

5390  Cherokee  Avenue 
Alexandria,  VA  22314 

1  Atlantic  Research  Corp. 

ATTN:  R.II.W.  Waesche 

7511  Wellington  Road 
Gainesville,  VA  22065 

1  AVCO  Everett  Research 
Laboratory  Division 
ATTN:  D.  Stickler 

2385  Revere  Beach  Parkway 
Everett,  MA  02149 


1  Battelle  Memorial  Institute 
Tactical  Technology  Center 
ATTN:  J.  Huggins 

505  King  Avenue 
Columbus,  OH  43201 

1  Cohen  Professional  Services 
ATTN:  N.S.  Cohen 

141  Channing  Street 
Redlands,  CA  92373 

1  Exxon  Research  &  Eng.  Co. 
ATTN:  A.  Dean 

Route  22E 

Annandale,  NJ  08801 

1  Ford  Aerospace  and 
Communications  Corp. 
DIVAD  Division 
Div.  Hq.,  Irvine 
ATTN:  D.  Williams 

Main  Street  &  Ford  Road 
Newport  Beach,  CA  92663 

1  General  Applied  Science 
Laboratories,  Inc. 

77  Raynor  Avenue 
Ronkonkama,  NY  11779-6649 

1  .  General  Electric  Armament 
&  Electrical  Systems 
ATTN:  MJ.  Bulman 

Lakeside  Avenue 
Burlington,  VT  05401 

1  General  Electric  Ordnance 

Systems 

ATTN:  J.  Mandzy 

100  Plastics  Avenue 
Pittsfield,  MA  01203 

2  General  Motors  Rsch  Labs 

Physics  Department 
ATTN:  T.  Sloan 

R.  Teets 
Warren,  MI  48090 

2  Hercules,  Inc. 

Allegheny  Ballistics  Lab. 
ATTN:  W.B.  Walkup 

E.A.  Yount 
P.O.  Box  210 

Rocket  Center,  WV  26726 
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1  Honeywell,  Inc. 

Government  and  Aerospace 
Products 

ATTN:  D.E.  Broden/ 

MS  MN50-2000 
600  2nd  Street  NE 
Hopkins,  MN  55343 

1  Honeywell,  Inc. 

ATTN:  R.E.  Tompkins 

MN38-3300 

10400  Yellow  Circle  Drive 
Minnetonka,  MN  55343 

1  IBM  Corporation 
ATTN:  A.C.  Tam 

Research  Division 
5600  Cottle  Road 
San  Jose,  CA  95193 

1  IIT  Research  Institute 

ATTN:  R.F.  Remaly 

10  West  35th  Street 
Chicago,  IL  60616 

2  Director 
Lawrence  Livermore 

National  Laboratory 
ATTN:  C.  Westbrook 

M.  Costantino 
P.O.  Box  808 
Livermore,  CA  94550 

1  Lockheed  Missiles  &  Space  Co. 
ATTN:  George  Lo 

3251  Hanover  Street 
Dept.  52-35/B204/2 
Palo  Alto,  CA  94304 

1  Los  Alamos  National  Lab 
ATTN:  B.  Nichols 

T7,  MS-B284 
P.O.  Box  1663 
Los  Alamos,  NM  87545 

1  National  Science  Foundation 
ATTN:  A.B.  Harvey 

Washington,  DC  20550 

1  Olin  Ordnance 

ATTN:  V.  McDonald,  Library 

P.O.  Box  222 

St.  Marks,  FL  32355-0222 


No.  of 

Copies  Organization 

1  Paul  Gough  Associates,  Inc. 

ATTN:  P.S.  Gough 

1048  South  Street 
Portsmouth,  NH  03801-5423 

2  Princeton  Combustion 

Research  Laboratories,  Inc. 
ATTN:  M.  Summerfield 

N.A.  Messina 
475  US  Highway  One 
Monmouth  Junction,  NJ  08852 

1  Hughes  Aircraft  Company 
ATTN:  T.E.  Ward 

8433  Fallbrook  Avenue 
Canoga  Park,  CA  91303 

1  Rockwell  International  Corp. 
Rocketdync  Division 
ATTN:  J.E.  Flanagan/HB02 

6633  Canoga  Avenue 
Canoga  Park,  CA  91304 

4  Sandia  National  Laboratories 
Division  8354 
ATTN:  R.  Cattolica 

S.  Johnston 
P.  Mattern 
D.  Stephenson 
Livermore,  CA  94550 

1  Science  Applications,  Inc. 
ATTN:  R.B.  Edelman 

23146  Cumorah  Crest 
Woodland  Hills,  CA  91364 

3  SRI  International 

ATTN:  G.  Smith 

D.  Crosley 
D.  Golden 

333  Ravcnswood  Avenue 
Menlo  Park,  CA  94025 

1  Stevens  Institute  of  Tech. 
Davidson  Laboratory 
ATTN:  R.  McAlevy,  III 

Hoboken,  NJ  07030 

1  Thiokol  Corporation 
Elkton  Division 
ATTN:  S.F.  Palopoli 

P.O.  Box  241 
Elkton,  MD  21921 
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1  Morton  Thiokol,  Inc. 

Huntsville  Division 
A'lTN:  J.  Deur 

Huntsville,  AL  35807-7501 

3  Thiokol  Corporation 
Wasatch  Division 
ATTN:  SJ.  Bennett 

AO.  Box  524 

Jrigham  City,  UT  84302 

1  United  Technologies  Research  Center 
ATTN:  A.C.  Eckbreth 

last  Hartford,  CT  06108 

3  United  Technologies  Corp. 

Chemical  Systems  Division 
ATTN:  R.S.  Brown 

T.D.  Myers  (2  copies) 

P.O.  Box  49028 

Sar  Jose,  CA  95161-9028 

1  Un  .  isal  Propulsion  Company 
ATTN:  HJ.  McSpadden 

Black  Canyon  Stage  1 
Box  1140 

Phoenix,  AZ  85029 

1  Veritay  Technology,  Inc. 

ATTN:  E.B.  Fisher 

4845  Millcrsport  Highway 
P.O.  Box  305 

East  Amherst,  NY  14051-0305 


1  University  of  California, 
Berkeley 

Chemistry  Deparment 
ATTN:  C.  Bradley  Moore 
211  Ixwis  Hall 
Berkeley,  CA  94720 

1  University  of  California, 

San  Diego 

ATTN:  F.A.  Williams 

AMES,  B010 
La  Jolla,  CA  92093 

2  University  of  California, 

Santa  Barbara 
Quantum  Institute 
ATTN:  K.  Schofield 

M.  Steinberg 

Santa  Barbara,  CA  93106 

1  University  of  Colorado  at 

Boulder 

Engineering  Center 
ATTN:  J.  Daily 

Campus  Box  427 
Boulder,  CO  80309-0427 

2  University  of  Southern 

California 
Dept,  of  Chemistry 
ATTN:  S.  Benson 

C.  Wittig 

Los  Angeles,  CA  90007 


1  Brigham  Young  University 
Dept,  of  Chemical  Engineering 
ATTN:  M.W.  Beckstead 

Provo,  UT  84058 


1  Case  Western  Reserve  Univ. 
Div.  of  Aerospace  Sciences 
ATrN:  J.  Tien 

Cleveland,  OH  44135 


California  Institute  of  Tech. 

Jet  Propulsion  Laboratory 
ATTN:  L  Strand/MS  512/102 

4800  Oak  Grove  Drive 
Pasadena,  CA  91009 


1  Cornell  University 

Department  of  Chemistry 
ATTN:  T.A.  Cool 

Baker  Laboratory 
Ithaca,  NY  14853 


1  California  Institute  of 
Technology 

ATTN:  F.E.C.  Culick/ 

MC  301-46 
204  Karman  Lab. 

Pasadena,  CA  91125 

1  University  of  California 
lots  Alamos  Scientific  Lab. 

P.O.  Box  1663,  Mail  Stop  B216 
Los  Alamos,  NM  87545 


1  University  of  Delaware 
ATTN:  T.  Brill 

Chemistry  Department 
Newark,  DE  19711 

1  University  of  Florida 
Dept,  of  Chemistry 
ATTN:  J.  Winefordner 

Gainesville,  FL  32611 
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3  Georgia  Institute  of 
Technology 
School  of  Aerospace 
Engineering 
ATTN:  E.  Price 

W.C.  Strahlc 
B.T.  Zinn 
Atlanta,  GA  30332 

1  University  of  Illinois 
Dept,  of  Mech.  Eng. 

ATTN:  H.  Krier 

144MEB,  1205  W.  Green  St. 
Urbana,  IL  61801 

1  Johns  Hopkins  University/APL 
Chemical  Propulsion 
Information  Agency 
ATTN:  T.W.  Christian 

Johns  Hopkins  Road 
Laurel,  MD  20707 

1  University  of  Michigan 
Gas  Dynamics  Lab 
Aerospace  Engineering  Bldg. 
ATTN:  G.M.  Facth 

Ann  Arbor,  MI  48109-2140 

1  University  of  Minnesota 
Dept,  of  Mechanical 
Engineering 
ATTN:  E.  Fletcher 

Minneapolis,  MN  55455 

3  Pennsylvania  State  University 
Applied  Research  Laboratory 
ATTN:  K.K.  Kuo 

K.  Palmer 
M.  Micci 

University  Park,  PA  16802 

1  Pennsylvania  State  University 
Dept,  of  Mechanical  Engineering 
ATTN:  V.  Yang 

University  Park,  PA  16802 

1  Polytechnic  Institute  of  NY 
Graduate  Center 
ATTN:  S.  Lederman 

Route  110 

Farmingdale,  NY  11735 


2  Princeton  University 
Forrestal  Campus  Library 
ATTN:  K.  Brezinsky 

I.  Glassman 
P.O.  Box  710 
Princeton,  NJ  08540 

1  Purdue  University 
School  of  Aeronautics 
and  Astronautics 
ATTN:  J.R.  Osborn 

Grissom  Hall 

West  Lafayette,  IN  47906 

1  Purdue  University 

Department  of  Chemistry 
ATTN:  E.  Grant 

West  Lafayette,  IN  47906 

2  Purdue  University 
School  of  Mechanical 

Engineering 

ATTN:  N.M.  Laurendeau 

S.N.B.  Murthy 
TSPC  Chaffee  Hall 
West  Lafayette,  IN  47906 

1  Rensselaer  Polytechnic  Inst. 
Dept,  of  Chemical  Engineering 
ATTN:  A.  Fontijn 

Troy,  NY  12181 

1  Stanford  University 
Dept,  of  Mechanical 
Engineering 
ATTN:  R.  Hanson 

Stanford,  CA  94305 

1  University  of  Texas 
Dept,  of  Chemistry 
ATTN:  W.  Gardiner 

Austin,  TX  78712 

1  University  of  Utah 

Dept,  of  Chemical  Engineering 
VTTN:  G.  Flandro 

^alt  Lake  City,  UT  84112 

1  Virginia  Polytechnic 
Institute  and 
State  University 
ATTN:  J.A.  Schetz 

Blacksburg,  VA  24061 
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1  Freedman  Associates 
ATTN:  E.  Freedman 

2411  Diana  Road 
Baltimore,  MD  21209-1525 


Intentionally  left  blank. 
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USER  EVALUATION  SHEET/CHANGE  OF  ADDRESS 


This  Laboratory  undertakes  a  continuing  effort  to  improve  the  quality  of  the  reports  it  publishes. 
Your  comments/answers  to  the  items/questions  below  will  aid  us  in  our  efforts. 

1 .  BRL  Report  Number _ BRL-TR-3165 _ Date  0f  Report  OCTOBER  1990 

2.  Date  Report  Reaped _ 

3.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related  project,  or  other  area  of  interest 

for  which  the  report  will  be  used.) _ 


4.  Specifically,  how  is  the  report  being  used?  (Information  source,  design  data,  procedure,  source 
of  ideas,  etc.) _ 


- £ - ;  ;  ■  ■  _ - - - - - 

-  \  .  p  «#■•.* ' 

5.  Has  tf§*  information  in' this  report  led  to  any  quantitative  savings  as  far  as  man-hours  or  dollars 
saved,  operating  costs  avoided,  or  efficiencies  achieved,  etc?  If  so,  please  elaborate. _ 


■  **\'7*r  “ 

Comments,  What  do  you  tHinlc  should  be  changed  to  improve  future  reports?  (Indicate 

•  changed  organization,  technical  content,  format,  etc.)  _ _ 

:  -  ;  '  ~~  - -  "" 

- 7 — ryg - - - — — -  -  — _ 

•  - 


Name 


CURRENT  Organization 

ADDRESS  _ 

Address 


City,  State,  Zip  Code. 

7.  If  indicating  a  Change  of  Address  or  Address  Correction,  please  provide  the  New  or  Coreect 
Address  in  Block  6  above  and  the  Old  or  Incorrect  address  below. 


Name 


OLD  Organization 

ADDRESS  _ 

Address 


City,  State,  Zip  Code 


(Remove  this  sheet,  fold  as  indicated,  staple  or  tape  closed,  and  mail.) 


•FOLD  HERE 


Department  of  the  Army 

Director 

U.S.  Army  Ballistic  Research  Laboratory 
ATTN:  SLCBR-DD-T 
Aberdeen  Proving  Ground,  MD  210b  -5066 
OFFICIAL  BUSINESS 


BUSINESS  REPLY  MAIL 

FIRST  CLASS  PERMIT  No  0001*,  APG,  MD 


POSTAGE  WILL  BE  PAID  BY  ADDRESSEE 


Director 

U.S.  Anny  Ballistic  Research  Laboratory 

ATTN:  SLCBR-DD-T 

Aberdeen  Proving  Ground,  MD  21005-9989 


NO  POSTAGE 
NECESSARY 
-  '  IF  WALED 
IN  THE 

UNITED  STATES 


-FOLD  HERE 


